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CONTACT PROCESSES 
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Cornell University 

Neuhauser [Probab. Theory Related Fields 91 (1992) 467-506] 
considered the two-type contact process and showed that on Z 2 co- 
existence is not possible if the death rates are equal and the particles 
use the same dispersal neighborhood. Here, we show that it is pos- 
sible for a species with a long-, but finite, range dispersal kernel to 
coexist with a superior competitor with nearest-neighbor dispersal in 
a model that includes deaths of blocks due to "forest fires." 

1. Introduction. In this paper we consider a two-type contact process in 
which the state of each site x £ Z 2 is = vacant, or 1 and 2 that indicate 
occupancy by two different species of trees. We formulate the evolution as 
follows: 

(i) Particles of type i die at rate 5i, give birth at rate 

(ii) A 1 born at x is sent to a y chosen according to a truncated power-law 
distribution pi(x,y) = c 1 3L { || J/ _ a .|| oo=1 } + c 2 l{ 1< \\ y - x \\ oc <M}\\y - x\\^ p , where 
Ci, C2 > and p < 3. 

(iii) A 2 born at x is sent to a y chosen at random from J\[\ = {y: \\x — 
2/||oo = 1}, that is, for these y, p2(x,y) = 1/8. 

(iv) If y is already occupied, then the birth is suppressed. 

(v) For each x, death to {y : \\x — y\\oo < F/2} occurs at rate 5q. 

The rules above imply that each type of particle behaves individually like 
a contact process, but in addition there are deaths to blocks of size F x F 
due to forest fires. 

Neuhauser [7] studied the model with <5o = and in which both species 
have the same dispersal distribution. She proved that coexistence is impos- 
sible on Z 2 if both species have the same death rate and conjectured that 
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in general the superior competitor (the species with the higher reproductive 
ratio Pi/Si) would win the competition. In this paper, we show that it is 
possible for two species to coexist in a model with forest fires if the weaker 
competitor has larger dispersal range. 

This result is, in a sense, an application of ecology to mathematics since 
it is based on a phenomenon that is well known to occur in forests; see [8]. 
The strongest competitors, which ultimately take over an undisturbed patch 
of forest, are those that need the least amount of light to grow. However, 
other species that are better at colonizing gaps in the forest created by 
disturbances can coexist with them. 

Our new result is not a counterexample to Neuhauser's conjecture, nor 
should it cast doubt on it. The competitive exclusion principle, which has 
been proved in the setting of ordinary differential equation models by Levin 
[5], states that the number of coexisting species cannot exceed the number 
of resources they compete for. (To be precise, Levin's result is for generic 
parameter values. In Neuhauser's model if (5\ = an d 5\ = 82, coexistence 
is possible in d > 3, but in this case we do not really have two species.) In 
Neuhauser's situation the two species compete for space, so there should 
only be one winner, even if different dispersal distributions are allowed. In 
our setting newly disturbed space is a second resource, so a species adapted 
to that niche can coexist with a superior competitor. 

To produce our example we begin by fixing /?2, 82 and the ratio Pi/ 81 so 
that 

(1) p 2 /(l + ft/fr) > W), P1/S1 > A C (M), 

where A c (Ai) is the critical value of (3/8 for the contact process with neigh- 
borhood J\f\. 

Theorem 1. Suppose (1). We can choose 8\, 8q, F and M so that Vs 
and 2's coexist. 

The second condition in (1) implies that when C2 = particles of type 1 
can survive in the absence of 2's. To explain the first condition, note that l's 
die at rate 8± and are born at a site at rate < Pi so the set of sites occupied by 
l's is dominated by a product measure with density Pi/ (81 + Pi). In the limit 
as 81 — > 00 with Pi/ 81 fixed, the environments seen by the 2's at different 
birth attempts are independent in the dominating process. Thus in the limit 
as 8\ — > 00 the 2's dominate a contact process with birth rate /?2<5i/(/3i + 8\) 
and death rate 82- Thus it suffices to have P2/O- + P1/81) > #2A c (A/"i) f° r the 
2's to survive in the system with no fire. Once this is done it is easy to use 
known results for the contact process to show if 8q < c/F 2 and c is small, 
2's also survive in the presence of fires. 
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Proposition 1. Suppose (1). There are constants c and A so that for 
any F > 0, the 2's survive whenever 5\ > A and 5q <c/F 2 . 

The l's survive by migrating from existing patches to newly created ones. 
By choosing the fire size F appropriately, the probability that a square 
contains a cleared area which remains open (free of 2's) for a long period 
of time is high. We choose 5q to produce enough such gaps, but also to not 
burn them again before we are ready. 

Proposition 2. Suppose (1). Fix 5\ > A. If <5 = a/F 3 with a suffi- 
ciently small, F > Fo(a) and M is chosen appropriately, then the l's sur- 
vive. 

Propositions 1 and 2 are proved by block constructions, so together they 
imply the existence of a stationary distribution in which both types are 
present with positive density. The remainder of the paper is devoted to 
proving the two propositions. 

2. Preliminaries on contact processes. In this paper, unless stated oth- 
erwise, we denote a contact process on Z 2 with neighborhood A/i and initial 
configuration A by Af. We let i,Af be the corresponding process truncated 
to not allow births outside [— L,L] d , and suppose that A C [—L,L] d . Define 
t a = min{t : Af = 0} to be the extinction time for Af, and lt A the extinc- 
tion time for L-Af. We now state two results for the contact process that 
will be used later on. First a survival condition, due to Bezuidenhout and 
Grimmett [1], which we take from Theorems 2.12 and 2.23 of [6]. 

Theorem 2. At survives if and only if it satisfies the following condi- 
tion: 

For every e > there are choices of n, L, T so that 

P{L+2nA [ T+i ]d D x + [-n,n] d for some x G [0,L) d } > 1 - e 

and 

P{ L +2nA+i n]i + [-n,n] d for some < t < T, 

and for some x G {L + n} x [0, L) d ~ 1 } > 1 — e. 

Note that in contrast to the usual block construction, one of the new 
copies of the completely filled square is on the side of the space-time box. 

The next theorem is a shape theorem proved by Durrett and Griffeath 
[3]. For the contact process they did this only for large enough birth rate. 
Given the results described in Chapter 2 of [6], it is routine to extend this 
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to all supercritical contact processes. Readers who do not want to take this 
leap of faith can replace A c (A/i) by the critical value for the one-dimensional 
nearest-neighbor contact process. 

Let rjt be a supercritical contact process, which includes the Richardson 
model (contact process with death rate 0) as a special case. Let v be the 
upper invariant measure of rjt, which is the limit starting from all l's, and 
let rjf be the stationary process starting from v. Let r = min{t : 77° = 0} and 
t(x) = min{t : x £ rjt}. The set of sites hit by time t is 

H t = {y G R d :3x £ Z d with Wx-yW^ < 1/2 and t(x) < t}. 

The region where rjf and rj\ are coupled is 

K t = {y eR d :3x£Z d with ||x-y||oo < 1/2 and rg(x) = Vt( x )}- 



Theorem 3. Let rjt be either a supercritical contact process or a Richard- 
son process (contact process with death rate 0) starting from 0. There is a 
convex set S C M. d such that for any e > 0, 

{l-e)S Ct- 1 H t C{l+e)S 

and 

(1 - e)S C r\H t n K t ) C (1 + e)S 
for t > To(e) a.s. on {r = 00}, where Tq(e) < 00 is a random time. 

3. Proof of Proposition 1. We begin by constructing our process from a 
collection of Poisson processes. For x,yeZ 2 and i = 1, 2, let {T% x ' v : n > 1}, 
{U^ x : n > 1} and {V® : n > 1} be the arrival times of Poisson processes with 
rates f3iPi(x,y), 5i and 6q. At times U^ x we put a <5j at rr to kill the particle 
at x if it is of type i. At times T^ x ' y we draw an arrow of type i from x to 
y to indicate that there is a birth from x to y if x is occupied by type i and 
y is empty. At time V x we put a 5 at all y such that ||y — x||oo < F/2 to kill 
all the particles inside the square. Our process has finite range so it follows 
from a result of Harris [4] that this gives a process well defined for all time. 

We first show that the 2's survive when there are no forest fires, that is, 
5q = 0. This is done by showing that the set of sites occupied by 2's can 
outcompete a process that is stochastically larger than the sites occupied 
by l's. Let £| be the set of sites occupied by i's when 5o = 0. Let Q be the 
set of sites occupied by i's when Sq = and each x € Z 2 not in state 2 flips 
from to 1 at rate (3\ and 1 to at rate 5\ independent of all other sites. 

Lemma 1. Suppose Co = Co- Using < for stochastically smaller than, we 
have & < C7- 
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Proof. Since the death rate of l's in ^ is equal to that of Q while the 
birth rate of l's is smaller, this follows from Theorem 1.5 in Chapter 3 of 
[6]. □ 

Lemma 2. Suppose 02, 82 and the ratio r = f3\/8\ are given. Let 7 = 
fo^i/iPi + £l) an d let A = Co- Writing a second superscript to indicate the 
value of 02, given e, 9 > there is a A such that there is a coupling with 

P{lA^- £ < L ( 2 < L A^ +£ forO<t<T)>l-9 

whenever 8\>A and A C 1? n [-L,L] 2 . 

Proof. In the graphical representation, each site x receives a type-2 
arrow at rate 02- The arrivals in [0, T] that touch some x G [—L,L] 2 are a 
Poisson point process. Let To = < T\ < T2 < ■ ■ ■ be the arrival times and let 
N(t) be the number of arrivals by time t. If N is large, P{N(t) > N) < 9/3. 
Having chosen this value of N, if a is small, 

Since rjt is a Markov process, the effect of a type-2 arrow to y at time t + h 
depends only on the state of the process at time t. The process that flips 
between 1 and at rates 8\ and (3\ is a two-state Markov process, so if <5i 
is large and /?i/<5i is fixed, the total variation distance between the process 
starting from 1 or and product measure is smaller than 9/3N at all times 
t > a. Since the process of 2-arrows is independent of the 1 flipping, and the 
latter has a stationary distribution (3\/{(3i +<5i), combining the last three 
estimates proves the result. □ 

Note that this implies ft =^ for all finite sets A and t < 00, where 
denotes weak convergence. Also by construction one has £f > Q f° r an 
t > 0. Therefore once we show that C,^ survives the forest fires, the proof is 
complete. But this is easy if we limit the frequency of our fires so that they 
are but small perturbations on the graphical representation. 

Proof of Proposition 1. For any chosen e, the idea is to choose Sq 
small enough so that the chance of a forest fire interfering with a given 
space-time block has probability < e and hence can be neglected. 

Let At be a contact process with birth rate 7' = fySi/iPi + Si) — e' and 
death rate 82- Since + ^1) > $AcCA/i)j we can pick e' small enough 

to have 7' > #2A c (-A/i). Then by Theorem 2, there are choices of n, L and T 
so that in the absence of fires 

P{ L+2nAj^™{ n ^ Dil [— n, n] 2 for some x G [0, L) 2 } > 1 — e 
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and 

P{L+2n^[+i^ ^ x + [ — n ' n ] 2 f° r some <t <T, 

and for some x & {L + n} x [0, L)} > 1 — e. 

Since a fire eliminates all the points within [x — F/2, x + F/2] 2 , our space- 
time block \—L — 2n, L + 2re] 2 x [0, T + 1] is affected only if the center of a fire 
is inside [-F/2 - L - 2n, F/2 + L + 2n] 2 . Let r = (F + 2L + 4re + 1) 2 (T + 1). 
The probability our space-time block is unaffected by fires is bounded below 
by 

1- f T 5 e- 5ot dt = e- 5oT > 1 -e 
Jo 

if ^ C /F 2 , where c > is some constant. Lemma 2 states that L+2nCt'^ " — 
L+2n^t U with probability > 1 — e whenever <5i > A. Combining the es- 
timates above we get 

P{L+2nCT [ +i' n]2 + \-n,nf for some x G [0,L) 2 } > (1 - e) 3 

and 

P{L+2n(?+i n ' n ^ 3i| [— n, n} 2 for some <t<T, 

and for some x G {L + re} x [0, L)} > (1 — e) 3 . 

Setting e small enough, we have recreated the block events of Theorem 1. 
For reasons that are explained in detail in Liggett's [6] book, £ 2 dominates a 
supercritical /s-dependent oriented percolation and therefore has a positive 
probability of survival. This completes the proof since by construction £f > 
Q 2 for all time. □ 

4. Proof of Proposition 2. 

4.1. Size of a fire box. We first determine an appropriate size F for our 
fires. Recall that the l's survive by colonizing the gaps created by forest fires. 
So we want a fire to create enough space for the l's to migrate from afar, 
and survive within the area for a long period of time. To accommodate these 
needs, we define three squares centered at 0, labeled B\, E>2 and B3 with 
sides of lengths 2L, 2r\L and 2r<iL = 2F, where 1 < r\ <V2- The innermost 
square B\ provides the space for l's to land, the square B2 allows them 
room to grow and the perimeter area B% — B2 serves as a buffer against the 
2's. We want to show the following. 

Proposition 3. Suppose 02, $2 are fixed and c\fi\ > 6i\ c (Afi). Here 
c\ is the coefficient for short-range births in the dispersal function pi(x,y), 
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and A c (A/"i) is the critical value of (5/5 for survival of the contact process 
with neighborhood M\ ■ Then there are choices of r\ , r<± and T so that in the 
absence of interference of other forest fires the following events happen with 
probability > 1 — 9 for large L: 

(i) If B\ receives a type-1 particle whose descendants survive within B2 
for a period of T, then they will survive within B2 for a period of 7T/2. 

(ii) If B3 does not contain any 2's at time 0, 2's do not reach B2 by time 
7T/2. 

In the first lemma we prove (i). To grow the set of l's we only use the 
short-range edges, and hence have a version of the process A x . Recall that 
Ti lt is the extinction time of Ti lA x . Let Hf and Kf , as in Theorem 3, be the 
set of sites hit and coupled for A x . For simplicity let Cf = Hf n Kf C Hf . 
Let Dji(y) be a square of side 2R centered at y and let Dr = Dr(0). Since 
(3 and 5 are fixed, there are i?i 5 2 > Ri,i > such that Dr 11 C S C Dr 12 , 
where S is the limiting shape of A® . Let 

(2) Gf = {(1 - e)D RlA c t-\C? - x),t-\Hf - x) C (1 + s)D Rh2 }. 

Lemma 3. Let At be a supercritical contact process with birth rate (5 and 
death rate 5. Let T = 8L/7R± : i and let x 6 [— L,L] 2 . For any 8 > 0, we can 
choose n so that for large L 

P{A x t C B 2 for t < 7T/2\ riL T x > T) > 1 - 6 

and 

P(r lL r x > 7T/2\ riL r x >T)>l — 9. 

Proof. Let e = 1/2 in (2). Theorem 2 implies that there is a t* > 
such that P(G? for all t > t*\r x = oo) > 1 - 6. For L > t*R 1A /A, we have on 

e = r\t>t* Gf, 

x s „ i r o r o rl2 r t rl2 



CI L/RlA ^x + [-2L,2LYD[-L,L] 



H x l/Ri i Ctx+ [-K, Kf where K = ±R^ 2 (4L/R 
Letting r\ = QRi,2/Ri,i + 1 we have on E, 

P([-L,L] 2 C C X L/Rl i C H X L/Ri i C D riL \r x = oo) > 1 
Now {G x t/2 ,t x = 00} C {G x t/2 ,t x > T}, so 

P(t x = 00) 
P(G x t/2 \t x >T)> P(G x T/2 \r x = oo) p l 
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The last fraction tends to 1 as T ^ oo. By the definitions of G x and r±, 

{r x >T}n G x T/2 = { riL r* >T}n G x 7T/2 . 
Since P( ri LT x >T)< P(t x > T) we have (recall 7T/2 = 4L/i? M ) 
P(G x 7T/2 \r x >T)< P(G x 7T/2 \ riL T x > T) 

<P( ri LT x >7T/2\ riL T x >T) 
and the desired conclusions follow. □ 

To prepare for the proof of Lemma 5 we observe that P(G 7T ^ 2 \t x = oo) — > 
1 as T — > oo so working backward through the last few calculations shows 
that for large L 

(3) P( riL T x >7T/2)>(l-0)P(r x = oo). 

The next lemma proves (ii) in Proposition 3. 

Lemma 4. If r 2 is chosen large enough, then, for large L, Aq n £3 = 
implies 

P(A t HB 2 = for all t < 7T/2) > 1 - 0. 

Proof. Fix e = 1/2. Let Rf be the Richardson process with neighbor- 
hood M\ and growth rate f3 2 - Let S be the convex set that is the limit in 
Theorem 3. Suppose R 2f2 > P%\ > are two numbers such that D R2 1 C S C 
Dr 2 2 . Let 

G t = {(1 - e)D R2A C r l B° t C (1 + e)D R2 2 }. 

There is a constant t* such that P(Gt for all i > t*) > 1 — 5. 
Suppose £2 > *i > Then 

P(G tl ) > P(G tl ,t 2 - 1 < C (1 + e)D R2 , 2 ) > 1 - 6. 

Since Rt is a Markov process, 

5>P(G tl ,fe-X^( 1 + £ ) i5 B 8 , a ) 
>P{t 2 ^RlzT DR2 ^(l + e)D R2a )P(G tl ). 
Since P(Gt 1 ) >l — 5, taking t 2 = t\ + 1 gives 

P((< + t!)- 1 ^ 1 ^ 1 ^ 2 ' 1 £ (1 + £ )^ 2 , 2 ) < */(l " S). 

Let 6» = 5/(1 - 6). Let m > r x L and m > (1 - e)t*i? 2 ,i- We can find t x > t* 
so that m = (1 — e)tii?2,i and as a result 

P(P^7 2 m]2 C (3/2)(7T/2 + h)D R22 ) >l-9. 
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J-i,- Jo,- J\- 



FlG. 1. Boxes in the block construction. 

Let r 2 = (3/2)(7T/2 + t\)R 2:2 /L and U = D T2 L- Since Rt is self-dual and 
B 2 C [—m, rrt] 2 , 

P{Rit/2 r\B 2 = 0) = P(R?t/2 nu c = 0) 
= P(R% /2 cU) 

>l-9. □ 

4.2. Block construction. Like Proposition 1, Proposition 2 will also be 
shown by using a block construction. We will make use of the following nota- 
tion. Let C = {(ttj, rt) € Z 2 : m + rt is even} and ei = (1, 0). We have already 
defined F = r 2 L. Let a > 0, k = F a , W = F + F 1+a , and define 

J_ = [-kW,0} 2 , J + = (0,kW} 2 , J mi ± = 2mkWe 1 + J±. 

We also define 

B(i,j) = [iW, (i + 1)W] x \jW, (j + 1)W), i,j e Z, 
B 2 {x,y) = [x — riL,x + riL] x [y - r\L,y + T\L]. 

We say there is a source of Vs at B(i,j) at stage rt if there is a block 
B 2 (x,y) C B(i,j) that is empty of 2's and contains at least one type-1 par- 
ticle at all times between 2nT and 2(n + 1)T, where T will be chosen later. 
We say there is a source of Vs in J mj + at stage n if there is a source of l's 
in one of the 5(i,j)'s that make up J m ,+- 

The next lemma is the key to proving Proposition 2. 

Lemma 5. Suppose we have a source of Vs at Jo,+ at stage 0. k and 
W can be chosen so that with probability > 1 — 9, there are sources of Vs at 
J_i _ and J\ _ at stage 1. 
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We alternate between + and — for our sources to avoid possible interfer- 
ence. Once Lemma 5 is proven, we will choose the cutoff M = AkW. The 
events described in Lemma 5 will have a finite range of dependence, so a 
standard argument (see, e.g., [2]) implies that our process will dominate a 
supercritical oriented percolation on C, and Proposition 2 follows. 

Proof of Lemma 5. There are three steps in creating a source in an 
adjacent box: 

(a) fire clears an area in [0, T/2] and no fire touches the clearing in 
[T/2,4T]; 

(b) dispersal brings a 1 to the box B\ inside the clearing during [T/2,T]; 

(c) the immigrant survives to time 4T without leaving B2, and without 
2's coming into B2. 

We will show that the intersection of these events has probability > 77 > 
in a given W xW square in J__i _ and Ji ; _. Since the J's are a grid of k 2 
such squares, if k is large, success will have high probability. At this point 
the reader might worry about 77 becoming small as k gets large but this 
point will be addressed in the proof. 

The probability of the first desired event happening in one specified W x 
W area is at least 

r (W-F) 2 T/2 r oo 

/ 5 e- 5ot dt- / 8 e- Sot dt 

JO J(2F) 2 7T/2 

= (1 - e -5o(W-F) 2 T/2^ e -l46 F 2 T _ 

If we choose W = F + F 1+a , a > 0, then to get 

e -14S F*T> 1 _ e and 1 _ e -6o(W-F) 2 T/2> 1 _ e 

we need 

<5 o T 2 T<-log(l-0)/14 and 5 F 2+2a T/2 > - log(0). 

The first is guaranteed by our assumption 5o < c/T 3 for a suitable choice of 
c. Once we fix 6q the second will hold when L and hence F is large. 

Suppose there is a source in Jo,+- It contains at least one particle for all 
times between T/2 and T. Then the maximum L°° distance between our 
source and an area cleared by a nice fire in J\- or J_i 5 _ (see Figure 1) is 
4kW. Thus, the minimum rate at which our source spreads to the target B\ 
within the cleared area is u = C2Pi\4:kW\~ p L 2 . If we let k = F a , then there 
is a c < 00 so that the probability of spread to the target area between T/2 
and T is bounded below by 

1 - exp(-uT/2) > 1 - expt-C/giL 3 -^ 2 ^). 
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To check this note that W < 2F +a and recall F and T are multiples of L. 
If (1 + 2a)p < 3, then for large L the last quantity is > 1 — 9. 

Once a particle has landed on an empty area, the probability that its 
descendant will survive there until t = 4T is bounded below by 

P(nZ^ T/2 ^0)>(l-0)P(r* = oo) 

as shown in (3). 

Combining the last three estimates, we see that the probability of pro- 
ducing a source in a given W x W area in J\ _ of J_i _ is bounded below 
as L — > oo. The number of independent opportunities, k = F a , tends to oo 
so the probability of success tends to 1 and the proof is complete. □ 
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